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p i' Abstract 



We calculate the enhancements to the inclusive charmless decays of b quark, 
b — > sg, sqq, sgg, from gluon penguin diagrams induced by the charged and neutral 
Higgs bosons {H^,hP ,H^ , and A9) in the Two-Higgs-Doublet Model with flavor- 
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C^ ' changing couplings. Within the considered parameter space, the new contributions 

from charged Higgs boson are dominant. After including the new contributions, 
the branching ratio BR{b -^ sg) {q^ = 0) can be increased form ~ 0.2% in the 
standard model to 4.4% and 2.6% in the two-Higgs-doublet model for mH+ = 100 
and 200 GeV, respectively. The new contribution to the decay mode h — > sqq with 
q = {u, d, s) is, however, numerically small and peaked at the lower g^ region. The 
new contribution to 6 — > sgg can also be neglected. 
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Among radiative decays of B meson, b ^ sg is theoretically clean, phenomenologi- 
cally interesting and sensitive to new physics beyond the standard model (SM), for ex- 
ample, the two-Higgs-doublet model (2HDM) [l[], the minimal supersymmetric standard 
model(MSSM) g § and Technicolor models^. 

In the SM, BR{b -^ sg) ~ 0.2% for on-shell gluon and BR{B -^ Xnocharm) ~ l-2%§. 
According to the studies in ref . ||^ , an enhanced h ^>- sgis favored phenomeno logically since 
such enhancement is very helpful for example to decrease the averaged charm multiplicity 
Uc and the semileptonic branching ratio[0 and to increase the kaon yields^]. For the 
large BR{B -^ v'-^s) measured by CLEO0, one can also give a plausible interpretation 
from an enhanced 6 — ^ s(7[^ [l^. The contributions to the ratio b -^ sg in 2HDM without 



flavor changing(FC) couplings were calculated in refs. [|ri|, |12| , and the authors found that 
this ratio less than 0.7% for mn > 200 GeV and tan/3 ~ 5 in both Type I and II 2HDM 
[0]. Such a enhancement is not large enough to meet the requirement. The possibility of 



a large rate for b —>■ sg in supersymmetric models were also studied in ref. |TB[ . 

In this letter, we calculate the contributions to the inclusive charmless decays b -^ sg 
(g^ = 0), 6 ^ sqq and b -^ sgg from the gluon penguin diagrams induced by the exchange 
of charged and neutral Higgs bosons in the so-called Model III: the two-Higgs-doublet 



model with PC couplings [p!4|, |T^. We found that the branching ratio BR{b —^ sg) can be 



increased from ~ 0.2% in the SM to 2 — 4% level in the Model III. So large enhancement 



is still consistent with the CLEO limit: BR{b -^ sg*) < 6.8% at 90%C.L.[|T6l, and will be 
very helpful to resolve the experimental puzzles^]. 

In the 2HDM, the tree level flavor changing scalar currents(FCSC's) are absent if 
one introduces an ad hoc discrete symmetry to constrain the 2HDM scalar potential and 
Yukawa Lagrangian. Lets consider a Yukawa Lagrangian of the form|T^ 

where (pi {i = 1, 2) are the two Higgs doublets, 0i,2 = '^T20i 2) Qi,L with i = (1, 2, 3) are the 
left-handed quarks, Uj^r and Djr are the right-handed up- and down-type quarks, while 
?7jj- and ^jj {i,j = 1,2,3 are family index ) are generally the nondiagonal matrices of 



the Yukawa coupling. By imposing the discrete symmetry (0i 
—Di, Ui -^ l^Ui) one obtains the so called model I and II. 
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In this letter, we will consider the third type of 2HDM: the so called Model III |T^ 
discrete symmetry is imposed and both up- and down-type quarks then have FC couplings 
with 01 and 02. In Model III, there are five physical Higgs bosons: the charged scalar 
H^, the neutral CP even scalars H^ and /i° and the CP odd pseudoscalar A^. After the 
rotation that diagonalizes the mass matrix of the quark fields, the Yukawa Lagrangian of 



quarks are the form [15 



^¥^ = v7,QLLhUj,R + r]EQi,L(j)lDj,R + &Qi,L(t)2Uj,R + ^FjQi,L<l)2Dj,R + h.C. 



(2) 



where 77^' = rriiSij/v correspond to the diagonal mass matrices of quarks and v ~ 246GeV^ 
is the vacuum expectation value of 0i, while the neutral and charged FC couplings will 



be [15, 17 
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where Vckm is the Cabibbo-Kabayashi-Maskawa mixing matrix[18[, and 



CU,D 



y/mffTj 
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In this letter we assume that Xij are real because we here do not consider the possible 
effects of CP violation induced by the phase of Ay. 



As pointed in ref. [15[ , the experimental data of K^ — K^ and B^ — B^ mixing processes 
put severe constraint on the FC couplings involving the first generation of quarks|T^: 
(A„j, Arfj) <^ 1 for i,j = (1,2,3). We here will enforce the same constraint: Xui^dj = for 
i,j = (1, 2, 3). And we also assume that Aj^ = 1, for i,j = (2, 3). 

Direct searches for Higgs bosons in 2HDM at LEP II[|19| place the following mass imits: 
Mh+ > 56GeV, M^o > 77GeV, Mao > 78GeV. From the CLEO data of BR{B -^ X,7), 



some constraint on M//+ in Model III can also be derived [|2y, ^ . If one uses new CLEO 
resuh [H of B{B -^ X,7) < 4.77 x 10"^ at the 3a level, the constraint Mh+ > 400 GeV 



can be read off directly from the Fig. 2 of ref.pOf. According to the studies in ref.[0], the 
existence of a charged Higgs boson with M^+ ~ 200 GeV is still allowed. In this letter, 
therefore, we will consider the mass range of 100 GeV to 800 GeV for all Higgs bosons in 
Model III. 

As for additional constraint on the Higgs boson masses from some other processes 
[ p^ , |2l[] , they are not good enough to compete with the F^ — F^ mixing processes. 

In the SM, the magnetic-penguin induced bsg coupling leads to b —>■ sg* transitions 
where g* could be light-like ( g^ = 0, on-shell gluon), or time-like (g^ > 0, off-shell gluon). 
Under the spectator approximation, there are basically three types of subprocesses: b — > 
sg, b -^ sqq for q = u,d, s, and b -^ sgg. It is straightforward to find the effective bsg 
coupling by explicit calculation, or by making appropriate changes to the sd'y vertex of 
ref . 1^ , 



with 



V,iq' 



r.(g^ 
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q Qpiu-qpiQy)! 



Fi^{q')L + F^{q')R 



+ta^,q'' [msFi{q^)L + m,F^{q^)R] } , (6) 

where p (ps) is the b (s) quark momentum, q = p — ps is the gluon momentum, F^ ' 
and F2 ' are the electric and magnetic form factors, gs ((72) is the QCD (electroweak) 
coupling constant, L,R = (1 =1= 75)/2 are the chirality projection operators and T" with 
a = 1, ■ ■ ■ , 8 are the SU{3)c generators. 

If the heavy top quark is the internal quark in the penguin diagram, the terms pro- 
portional to ml/M^r, ml/My^r, q^ /^^'i can be neglected and we then have[^, ^ 

^i^(O) = Y.^^h{x,), F/^(0) = 0, (7) 

i 

F,\0) = F,^{0)=j:nMx,), (8) 



with 



fi{xt) 



i2{i-xty 



18xt - 29x,2 + lOxj* + xf-(8- 32xt + 18x?) Inx* 



(9) 



f2{xt) = -rr, ^7 2 + 3xt-Qx'f + x^t+6xt\nxt , (10) 

4(1 -XtP L J 

where Xt = m^/Myy, vi = V*.Vib for i = u,c, t. 

For decay b ^ sg, q^ax ~ 20Gey^ and the assumption q^ <^ mf which would justify 
the replacement of the form factors with their values at g^ = is valid only for top quark, 
but not for the light u and c quarks ^, ^. The correct fi{xj) for j = u,c have been 
given for example in ref.|0, and we also calculated the f2{xj) for j = (m, c) by using the 
same technique as refs.|^, |25[| , 

,r 2^ 10 2, 2 2(1 + 2z,) , , ,^^, 

f2ixj,q^) = -Xj{\nxj + 2g{zj)), (12) 

where Zj = g^/(4m|), and the explicit expression of g{zj) can be found for example in 
ref. 0. For g^ > 4m|, the internal u or c quark are on mass-shell, an absorptive part 
therefore appears for fi{xj) and f2{xj), respectively. 

In the Model III, an effective bsg coupling can also be induced by the gluon penguin 
diagram via exchanges of charged and neutral Higgs bosons {H^, h^,A^, H^) as depicted 
in Fig.l. We will evaluate the penguin diagrams involving H^ and then extend the 
calculation to the neutral Higgs bosons. We will follow the same procedure as that in the 
SM. We will use dimensional regularization to regulate all the ultraviolet divergence in 



the virtual loop corrections and adopt the MS renormalization scheme. The necessary 
Feynman rules can be obtained from the Lagrangian in Eq.(|D. 

The Feynman diagrams with light internal u and d quarks do not contribute since we 
have assumed that Xui^j = for z, j = (1, 2, 3). By explicit analytical calculations, we can 
extract out the form factors F^ ' and Fg ' which describe the new contributions from 
the neutral and charged Higgs bosons. 

For the case of charged Higgs boson H^, the effective bsg vertex is of the form, 

,2 

2_ 
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r.(g^) = ^,^MPs)T^{{q'9,.'q,qu)Y[F,'L + F,^R 
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(14) 
(15) 
(16) 
(17) 



with 



Ci{xi) 
Cii{x,i) 
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logXi 
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36(1 -x,)3 ' 6(l-x,)4^°^^'' 
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(18) 
(19) 



^^ = y(Vb.Me")l(^c.Me")^,, 5. = ^(eVc.M);(^^V^. 
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C. = -!L{VcKMe)l{eVc,M)„, D, = '"-.{fVc,,,)l{eV, 



'CKMj.^i^O) 



where ^'^'■^ have been given in Eq.(|^), Xi = mf/Ml+ for i = c,t. The expressions of 
Ci(xj) and Cii(xj) function are exact for the case of heavy internal top quark, and correct 
approximately ( the error is less than 5%) for the case of internal c quark. 

By the same procedure, we get the effective bsg vertex induced by the neutral Higgs 
bosons, 



where the summation over all three kinds of neutral Higgs bosons is understood. 
For the CP even scalar Higgs boson h^, we have 

1 



(21) 
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^ [3Cn(x.)-Ci(a;,)] 



nh" 



^Jm^rribmiXisXib cos a — y/msmbinbXbs sin a cos a 
1 



^Ml. ,=,, 



^ [-mb(3Cii(xi) - Ci{xi)) - 2miCi{xi)\ 



^/rn^rnbniiXisXib cos^ a — y/m^mbnibXbs sin a cos a 



(22) 



(23) 



where Xi = m^/M^o- 

For the CP even scalar Higgs boson H^, we have 



-ff" i=s,b 



:^H0 



y/m^mbTniXisXib sin a + y/msmbmbXbs sin a cos a 
TnTT Yl [-"^b(3Cii(xj) - Ci{xi)) - 2miCi{xi)] 



^/rn^rnbrriiXisXib sin^ a + y/m^nib'i^hXbs sin a cos a 

where Xj = mj/Mfjo- 

For the CP odd pseudoscalar Higgs boson A^, we have 



(24) 
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Y^ [3Cu{xi) - Ci{xi)] y/m~mbmiXisX 



ibi 
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^ [-nib (3Cii(xj) - Ci(xi)] 



-2mjCi(si)] y/msm^niiXisXib, 



(26) 



(27) 



where Xj = mf/M\o, and the functions Ci{xi) and Cii(xj) in Eqs.( P^p7D have been given 
inEqs.([gjT|). 

In the numerical calculations, we fix the following parameters and use them as the 
standard input (SIP) ^^■. Mw = SOAlGeV, a^^ = 1/129, sin^ ^vk = 0.23, m„ = 5 



MeV, nid = 9 MeV, m^ = 1.4GeV, m, = 0.13 GeV, mb = 4.4GeV, mt = 170GeV, 
Ag_ = 0.225GeV, A = 0.84, A = 0.22, p = 0, and r/ = 0.36. For the definitions and values 
of these input parameters, one can see refs. p6t ^% . For the running of a;s(/i), we use the 



two-loop formulae as given in ref. [^ and take /i = nib- 
Using the SIP and assuming muiggs = 3006*61^ and g^ = (1 — \9)GeV'^, we find 
numerically that: 



;ii,H| 



In the SM, \Fl\ ^ 0.25, iFg"^'"! ^ 0.008; For all neutral Higgs bosons, Fi and F2 
form factors are less than 10"'^ in magnitude and therefore can be neglected safely; 



\L.R 



For charged Higgs boson, F^ ' are very small and therefore can be neglected, but 



i^ ~ 1-2 is larger than its SM counterpart by about two orders of magnitude. While 
the F^ ^ 0.009 is comparable with F^ in the SM. 

In the following, we will only consider the new effects of Fg ' form factors induced 
by charged Higgs boson. The q dependence of Fi and F2 in the SM and the F2 ' in the 
Model III comes from the penguin diagrams with c internal quarks. But this dependence 
is very weak and therefore can be neglected safely. 

Now we turn to calculate the inclusive charmless decay b ^ sg for both on-shell and 
off-shell (g^ > 0) gluon, in order to check the possible enhancements to the decay widths 
and branching ratios induced by the charged Higgs boson in Model III. 

For the decay b —* sg with an on-shell gluon as depicted in Fig. 2a, it is easy to derive 
the decay rate 



T{b ^ sg)'^' = ^^ "^ Fo \F2\' (2^ 

TT 



in the SM, and 



r(b^ss)- = ?^^^Tj\F,r-+'4m') (29) 

TT V ^b J 

in the Model III. Here Fq = G^mf/(1927r3), F^'^ = F2 + F^'^, the SM form factor F2 
has been given in Eq.(RD. For M^+ = 200 GeV, one finds 



F(6 ^ sgY^^ 
T{b-^ sg)' 



^^ = r . ^ "sM = 14.6 (30) 



and the mass dependence of Rg is shown in Fig. 3a 
The branching ratio BR{b — » sg) is of the form 



3a,(mfe) -0.105 



IF2P, m SM, 



(31) 
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where f{z) = 0.54 is the phase space factor in the semileptonic b-decay with z = rric/mh, 
k{z) = 0.88 is the QCD factor to the semileptonic b-decay, and BR{B —>■ XcCVg) = 
10.5%pl. 



Figs. (3a) and (3b) are the plots of Rg and the ratio BR{h -^ sg) versus Mh+ in all 
three types of 2HDM for Mh+ = (100 - 800) GeV. For theoretical predictions in Model 



I and Model II, we directly used the formulae given in ref.|TT|. For Mh+ = 200GeV, the 
branching ratio BR{b -^ sg) will be increased from ~ 0.2% in the SM to ~ 2.6% in the 
Model III, as shown by the upper solid curve in Fig. 3b. The short-dashed and dotted 
curve in figs. (3a) and (3b) shows the ratios in Model I and Model II, respectively. It is 
easy to see that the enhancements to the ratios in both Model I and II are much smaller 
than that in Model III. 

For g^ > 0, the gluon is virtual and time-like, therefore before it can fragment into 
real hadrons, it first disintegrates into real partons such as on-shell quark pair qq for 
q = {u, d, s) and gluon pair gg, as shown in Fig.2. The process b ^ uu and b ^ dd 
can be treated on the same footing, while one should take into account the identical 
particle effects for the decay b —>■ sss. Here, as an illustration, we will evaluate the new 
contributions to the branching ratio BR{b -^ sqq) for g = u, d, s in the Model III. We use 
the same method as ref . |p^4|| . 

In order to illustrate the effects of charge Higgs boson on the differential decay rate, 
we draw the Fig. 4 for the typical decay mode b — > suu. 
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1 - y)Re {F^YF^ , (32) 



where y = q^ /rn^ ^ 1 — 2Es/m,i,. If one treats the Eg as the kaon energy, then the 
differential rate in Eq.(^) can be regarded as the " Kaon-energy" spectrum. In Fig. 4, 
the cc threshold cusp is clearly exhibited, where the short-dashed (solid) curve shows 
the differential rate in the SM (Model III). It is easy to see that the new contribution is 
peaked at the lower g^ region. 



The decay width r(6 — * sqq) with q = {u, d, s) is of the form 
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where x^m = ymin = 4mg/m^. Contrary to the case of decay b —* sg with on-shell gluon, 
the new contribution from charged Higgs boson in Model III tend to decrease the decay 
width r(6 — > sqq) slightly by about five percent with respect to the SM prediction. The 
reason for this behaviour is simple. We know that only the new magnetic F2 form factor 
contribute effectively to the decay processes under study. For the decay b ^ sg with 
on-shell gluon, the magnetic form factor F2 dominate the total contribution. But for 
the case of g^ > 0, the SM form factor F^ control the decay processes, the new F2 can 
contribute only through interference with the SM F(". The large infrared(IR) logarithms 
in Eq.(|33|) will be canceled if one make a complete 0{as) treatment of the second term in 
Eg. (|53D . Since we are not giving a full 0{as) QCD analysis, we can not include this term 



consistently. We therefore simply drop it from further discussions, as done in ref. [^ . 

The decay mode b — > sgg via or not via g* —>■ gg has been studied in ref. |^^. It was 
found that the decay width r(6 — > sgg) depend on the form factor Fi only, 
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As discussed previously, the new contributions to the electric form factor Fi from charged 
and neutral Higgs bosons are very small and has been neglected. The decay b -^ sgg there- 
fore will be not affected effectively by gluon penguins induced by Higgs bosons appeared 
in Model HI. 

Collectively, the branching ratio BR{b —>■ sg*) (here b — > sg* is symbolic for processes 
oi b -^ sg (on-shell gluon), b -^ sqq, and b -^ sgg ) can be written in the form 
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As illustrated in Fig. 5, the branching ratio BR(h -^ sg*) = 5.5% — 1.6% for Mh+ = 
100 — 800 GeV in the Model III (solid curve) compared with 1.3% in the SM (short-dashed 
line). For light charged Higgs, the enhancement is significant: BR{b —>■ sg*) = 3.7% for 
Mh+ = 200 GeV in the Model III. The upper dot-dashed line in Fig.5 shows the CLEO 



upper limit: BR{b -* sg*) < 6.8% at 90%C.L.[||]. In ref.[|8[, the author studied the CP 
violation in radiative B decays in a type-Ill 2HDM with the fourth generation fermions, 
and found that the branching ratio BR{b -^ sg) can be increased to 10% level. 

To summarize, we have calculated, from the first principle, the new gluon penguin 
diagrams that contribute to the inclusive charmless decays of b quark in the Model III 
without inclusion of QCD corrections. We start from the evaluation of the new gluon pen- 
guin diagrams induced by the exchange of charged and neutral Higgs bosons {H^, h^, H^, 
and A^), derive out the Ri and R2 form factors which control the new contributions to 
the inclusive b quark decays under study, and finally calculate the relevant decay rates 
and branching ratios. We found that: 

(a) Among charged and neutral Higgs bosons, the CP even charged Higgs 
boson H^ dominate the new contribution. By using the assumed FC cou- 
plings, the contribution of the neutral scalar and pseudo-scalar is completely 
negligible. Therefore, both the value of the mixing angle a and masses 
{mhO,mHO, and rri^o) are irrelevant. 

(b) The new electric form factor Ri is also completely negligible. The new 
magnetic form factor i^, on the contrary, is much larger than its SM coun- 
terpart and thus contribute significantly to the decay processes in question. 
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The F-^ is comparable in size with the F2 in the SM and can also contribute 
effectively. 

(c) The charged Higgs enhancement to the branching ratio BR{b -^ sg) 
with g^ = can be as large as a factor of 25 (14.6) for Mh+ = lOOGeV" 
{200GeV) in the Model III. So large enhancement will be very helpful to 
generate a large ratio BR{b -^ sg) favored by some experimental data such as 
the deficit in charm counting, a 3a deficit in kaon counting as well as the well- 
known large B -^ fj'Xg branching ratio measured by the CLEO collaboration. 
The enhancement in the Model III can be much larger than that in both Model 
I and II. 

(d) For the subprocesses with a off-shell gluon (g^ > 0), the new contri- 
bution to the decay mode b — > sqq with q = {u, d, s) is numerically small and 
peaked at the lower q'^ region. And the new contribution to 6 ^ sgg can be 
neglected. 
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Figure Captions 

Fig.l: Self-energy and gluon penguin diagrams witli cliarged and neutral scalar or pseu- 
doscalar exchanges in the Model III. The charged and neutral Higgs boson prop- 
agators correspond to the up-type quarks Ui = {u, c, t) and the down-type quarks 
di = {d, s, b), respectively. 

Fig. 2: Feynman diagrams for the charmless decay of b quark in the Model III: (a) b ^ sg 
(g^ = 0, on-shell gluon); (b) b -^ sqq (g^ > 0) with q = u,d, s; (c) additional diagram 
for b — * sss;(d) b -^ sgg via g* -^ gg; (e) typical one-particle reducible diagram that 
lead to b ^ sgg; (f) typical one-particle irreducible diagram that lead to 6 — > sgg. 
The blobs in first five diagrams denote the effective bsg vertex. 

Fig. 3: Mass dependence of the ratios Rg and BR{b -^ sg) {q^ = 0) in three type of two- 
Higgs-doublet models for Mh+ = 100 — SOOGeV. In Fig. 3a and (3b), the dotted, 
short-dashed and solid curve shows the ratios in the Model I, II and III, respectively. 
The dot-dashed line in Fig. 3b is the SM prediction. 

Fig. 4: The g^ dependence of differential decay rates for b — > suu. The short-dashed and 
solid line shows the rate in the SM and Model III, respectively. 

Fig. 5: The mass dependence of branching ratio BR{b -^ sg*) for Mh+ = 100 — SOOGeV. 
The short-dashed and solid curve shows the ratio in the SM and the Model III, 
respectively. 
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